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AGN Reverberation Mapping 
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Abstract Reverberation mapping is now a well-established technique for investi¬ 
gating spatially-unresolved structures in the nuclei of distant galaxies with actively- 
accreting supermassive black holes. Structural parameters for the broad emission- 
line region, with angular sizes of microarcseconds, can be constrained through the 
substitution of time resolution for spatial resolution. Many reverberation experi¬ 
ments over the last 30 years have led to a practical understanding of the require¬ 
ments necessary for a successful program. With reverberation measurements now in 
hand for 60 active galaxies, and more on the horizon, we are able to directly con¬ 
strain black hole masses, derive scaling relationships that allow large numbers of 
black hole mass estimates throughout the observable Universe, and begin investi¬ 
gating the detailed geometry and kinematics of the broad line region. Reverberation 
mapping is therefore one of the few techniques available that will allow a deeper 
understanding of the physical mechanisms involved in AGN feeding and feedback 
at very small scales, as well as constraints on the growth and evolution of black 
holes across cosmic time. In this contribution, I will briefly review the background, 
implementation, and major results derived from this high angular resolution tech¬ 
nique. 


1 Introduction and Motivation 


In the 25 years that the Hubble Space Telescope has been in low Earth orbit, there 
have been leaps and bounds in our understanding of many astrophysical phenom¬ 
ena, not the least of which are supermassive black holes. Dedicated ground-based 
programs coupled with the exquisite spatial resolution afforded by HST have led 
to the now-common understanding that massive galaxies host supermassive black 
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holes in their cores (see the review by Ferrarese & Ford 2005). In general, studies 
have shown that more massive galaxies host more massive black holes, which is 
interpreted as a symbiosis between galaxies and black holes in which they grow to¬ 
gether and regulate each other’s growth. Recent studies have begun to find several 
galaxies that do not appear to follow this simple scaling relationship, however, and 
our picture of galaxy and black hole co-evolution throughout the Universe’s history 
is becoming more complicated (see the review by Kormendv & Holl2013l) . 

The active galactic nucleus, or AGN, phase is now understood to be a sporadic 
event in the life of a typical supermassive black hole, thought to be t riggered by a 
mer ger or secular process in the host galaxy (see the review by Hec kman & Best 
,20141 and references therein). During this phase, the black hole is accreting at a 


relatively high rate, and the accretion process is releasing large amounts of energy 
across the electromagnetic spectrum. The lifetime of a typical AGN event for a 
typical black hole is small compared to the age of the Universe, and is generally 
thought to be on the order of ~ 10 8 years (e.g., Kelly et akll2010l) . Cosmic downsiz¬ 
ing is observed in AGNs as it is in galaxies: the bright quasars we see in the early 
Universe are associated with massive black holes and large accretion rates, while to¬ 
day’s active black h oles tend to be fewer in number, less massive , and have smaller 
accretion rates (e.g., Ueda et al. 2003 . 2014 : Shankar et al. 2009t Kelly et al. 2010t 


Kalfountzou et al. 2014) . 

Unfortunately, the rarity of bright AGNs compared to the plethora of galaxies 
hosting quiescent black holes today leads to the situation in which we find ourselves, 
where AGNs are generally distant and difficult to study even with the pristine spatial 
resolution afforded by HST. The techniques that have been developed to constrain 
the masses of inactive black holes in nearby galaxies rely on spatially resolving the 
innermost parsecs of a galactic nucleus, so they are limited to distances < lOOMpc 
and therefore not applicable to most AGNs. And yet AGNs act as beacons that shine 
across the entire observable Universe, tempting us to try to understand black hole 


2014j) and perhaps 


and galaxy growth and evolution out to z = 7.1 (iMomiian et al 
beyond. 

Luckily, AGNs are not only bright, they are also highly variable. We can, in 
effect, substitute tim e resolution for spatial resolution in a tech nique known as re¬ 
verberation mapping ( Blandford & McKee! 1982 ; Petersonll 19931) to probe microarc- 
second scales in the nuclei of even the most distant active galaxies. 


2 Reverberation Mapping Primer 


Many independent studie s have led to the general picture of AGN structure that 
we understan d today (e.g., Antonuccill 1 993b Urrv & Padovanifl 995 : Petersonll 1997 : 
NetzerlEo 1 5l and references therein), as represented by the cartoon diagram in Fig- 
ureQ] In the center is the supermassive black hole, with a mass in the range of 
10 6 — 10 10 Mq, and its associated accretion disk. The jet (if the AGN has one) is per¬ 
pendicular to the accretion disk and highly collimated. On larger scales (~ 0.01 pc 
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Fig. 1 Cartoon diagram of the 
typically-assumed structure 
of an AGN. Depending on 
the orientation at which this 
stmcture is viewed, different 
spectral signatures will be 
seen. In particular, the BLR 
is only visible to observers 
with a relatively face-on view. 
If viewed from the side, the 
torus blocks the BLR and 
only narrow emission lines 
will be seen in the AGN spec - 
trum. From lUrrv & Padovanil 
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for typical Seyferts, approximately the extent of the inner Oort Cloud in our own 
Solar System) lies a region of photoionized gas that radiates line emission. The lo¬ 
cation of this gas deep within the potential well of the black hole results in line-of- 
sight gas velocities that are quite large, causing the emission lines to appear Doppler 
broadened in the AGN spectrum by a few 1000 km s“*. We imaginatively call this 
region of gas the broad line region (BL R). The outer edge of the BLR is most likely 


set by the dust su blimation radius (e.g.. lNetzer & Laonll993l : lNenkova et al.ll2008 


Goad et al.lEom as the inclusion of dust in the BLR gas will extinguish the line 


emission. The dusty gas that exists outside this radius is generally referred to as 
the “torus”, although the exact geometry of the region is not known. The dust torus 
causes the AGN system to have different spectral signatures when viewed at differ¬ 
ent orientations — a system that is close to face-on will have broad emission lines in 
its spectrum, while a system that is viewed edge-on will have the dust torus blocking 
the observer’s line of sight to the BLR, so no broad emission lines will be seen. On 
even larger scales (~Tens of pc for typical Seyferts), additional gas that is photoion¬ 
ized by the AGN system also exists, but the location of this gas on galactic scales 
results in line-of-sight velocities that are on the order of a few 100 k m s 1 . We see 
the signature of this gas as narrow emission lines in the AGN spectrum, hence the 
name attributed to this region of gas is the narrow line region (NLR). 

Reverberation, or “echo”, mapping measures the light travel time between dif¬ 
ferent regions in an AGN system. The continuum emission is expected to arise 
from the accretion disk, and in the ultraviolet, optical, and near-IR it is observed 
to vary on timescales of hours to days. The source of this variability is not yet un- 
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JD (-2,400,000) Time delay (days) 


Fig. 2 Continuum and H/3 light curves for the Seyfert galaxy NGC4151 (left), autocorrelation 
function of the continuum light curve (top right), and cross-correlation of Hj8 relative to the con¬ 
tinuum (bottom right). The emission-line light curve is visibly delayed from that of the continuum 
and smoothed in time, both signatures that are obvious to the eye in the light curves and the cross¬ 
correlation function, and indicative of the extended nature of the Hj8-emitting BLR gas. From 
Be ntz et all l2006l) . 


derstood (one possib le explanation is magnetic recombination in the accretion disk, 
e.g., Kawaguchi et all2000l) . but whatever the cause, the variations that are observed 
in the continuum emission are seen echoed at a later time in the fluxes of the broad 
emission lines (see Figure[2]for example light curves for the nearby Seyfert galaxy 
NGC4151). The time delay between the variations in the continuum and the echo 
of those variations in an emission line is simply the average light travel time from 
the accretion disk to the BLR. The accretion disk is generally assumed to be very 
compact, and so the time delay can be interpreted as the average radius of the BLR 
in the AGN system. By definition, reverberation mapping requires a line of sight 
that permits the observer to view the broad emission lines in the AGN spectrum, so 
it is only applicable to AGNs with a relatively face-on orientation. 

Within the BLR, different emission lines are observed to respond to continuum 
variations with different time delays, such that species with higher ionization poten¬ 
tials, like C IV A1549, respond with a shorter time delay than those with lower ion¬ 


ization potentials, like H/3 (e.g., Peterson & Wandel 2000t Kollatschnv et al. 2001 


Bentz et akll2010bl) . This behavior points to ionization stratification within the BLR 
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- more highly ionized line emission is radiated from a smaller radius within the 
BLR, while emission from more neutral gas occurs at larger radii, further from the 
central ionizing source. Photoionization modeling constraints agree with this inter¬ 
pretation: the photoionized gas in the BLR preferentially emits line emission wher¬ 
ever the temperature and density are most favorable for a specific at omic transition 
(e.g., Baldwin et all 1995 : Korista et~ail 1 997t Korista & GoadlEoQ4 ). 

With reverberation mapping, the radius we measure is the responsivity-weighted 
average radius, which may not exactly coincide with the emissivity-weighted av¬ 
erage radius of the gas for that particular line emission. But the general behavior 
that we observe is the same: in a single AGN, if the central luminosity increases, 
then the time delay we measure for any specific emission line be comes larger (e.g., 
Peterson et al. 2002 : Bentz et all 2007 : Kilerci Eser et al.l 2015b . It is important to 
keep in mind that we are not measuring a geometrical radius, such as an inner or 
outer boundary of the BLR, with emission-line reverberation mapping. 

Even for a relatively nearby AGN, the size of the region probed by reverberation 
mapping is quite compact and not resolvable with current imaging technology or 
that which is likely to be developed in the foreseeable future. For a typical Seyfert 
galaxy at a distance of 40 Mpc, the time delay expected for the Up A4861 emission 
line is ~ 10 light-days, which projects to an angular radius of ~ 50microarcsec. 

In practice, reverberation mapping relies on dense spectrophotometric monitor¬ 
ing of an AGN system over an extended period of time. The basic requirements for 
a successful monitoring program are the following: 

1. total campaign length at least three times the longest expected time delay to max¬ 
imize the probability th at v aribility of a large enough amplitude will occur during 


the campaign ( Home et al. . 20041) : 


2 . 


sampling cadence that is sufficiently dense to resolve the variability and expected 
time delays (e.g., Nyquist sampled in time); 

3. exposure times that yield a signal-to-noise > 50 in the continuum and substan¬ 
tially higher in the emission lines, where the amplitude of variability is generally 
only a few percent; 

spectral resolution that is sufficiently high to distinguish broad emission lines 
from each other and from overlapping or nearby narrow emission lines (generally 
R of a few thousand); 

flux calibration that is good to 2% or bette r from observation to observation 


4. 


5. 


throughout the campaign (IPeterson et al .. 2004); 


6. strong nerves and a healthy dose of good luck. 

Meeting all of these constraints with a ground-based telescope is observationally 
quite challenging. The weather especially can cause an otherwise well-planned re¬ 
verberation campaign to not live up to its potential or to fail outright. Added to this 
uncertainty is the fact that AGN variations are stochastic and not guaranteed to oc¬ 
cur during the course of a monitoring campaign for any particular AGN of interest 
(cf. the case o f Mrk 290 which wa s monitored in 20 07 and showed stro ng variations 
as reported by Penney et al.ll2Qlol but according to Bentz et akll2009b showed little 
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to no variation when it was monitored again in 2008). The typical warning for the 
stock market also applies here: past performance is no guarantee of future behavior. 

To date, s ucces sful reverberation campaigns have been carried out for ~ 60 dif¬ 
ferent AGNs ( Bentz & Katzl 2015 ). Constraints on readily-available resources have 
generally limited the size of the telescope used in a reverberation campaign to 
1.0 — 4.0-m class telescopes. Coupled with the need for high signal-to-noise ratio 
spectra in each visit, most of these 60 AGNs are apparently bright and reside within 
the nearby Universe (z < 0.1). Reverberation mapping is not fundamentally lim¬ 
ited to nearby objects, but the high luminosities necessary for high signal-to-noise 
spectra of z = 2 — 3 quasars directly translates to a long emission-line time delay, 
which is further stretched through cosmological dilation (cf. Kaspi et akl 2007 and 
their monitoring campaign length of 10 years). Ongoing efforts to multiplex rever¬ 
beration mapping with the Sloan Digital Sky Survey multi-object spectrograph may 
soon increase the sample size by a s ubstantia l fraction and push the median of the 


sample to somewhat larger redshifts (IShen et all 120151). but large statistical samples 


of reverberation results at multiple redshifts spanning the course of cosmic history 
are unlikely for the forseeable future. Luckily, the results we have in hand for the 
current reverberation sample are able to provide us with a foothold for investigating 
cosmological black hole growth. 


3 Reverberation Mapping Products 

Substantial progress over the last 10 years especially has led to several valuable 
and widely used reverberation-mapping products, including direct black hole mass 
measurements, black hole scaling relationships that can be used to quickly estimate 
large numbers of black hole masses, and detailed information on the geometry and 
kinematics of the BLR gas. As we describe in the following sections, the first two 
items are fairly well-developed at this time, while we are just now beginning to truly 
exploit the third. 


3.1 Black Hole Masses 


In order to directly constrain the mass of a black hole, a luminous tracer (usually gas 
or stars) must be used to probe the local gravitational potential. In the case of rever¬ 
beration mapping, the photoionized BLR gas is deep within the potential well of the 
black hole where its motion is expected to be dominated by gravity in the absence 
of strong radiation pressure. While the effect of radiation pressure is still debatable 
for lines such as C IV, these condi t ions are most likely to be satisfied for H/3 (e.g., 
Marconi et al. 2008 : Netzer 20091 Netzer & Marzianil 2010l) which is expected to 
arise from gas that is well-shielded from the central ionizing source. The Doppler- 
broadened width of the H/l emission line is therefore a measure of the line-of-sight 
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velocity of the gas within the BLR. And, as described in Section[2] the time delay 
in the Hj3 emission line is a measure of the radius of the BLR for that same line- 
emitting gas. Coupling these two measurements together through the virial theorem 
allows a direct constraint on the mass of the central black hole, modulo a scaling 
factor to account for the detailed geometry (including inclination) and kinematics 
(whether rotation, infall, or outflow) of the gas. 

Given the orientation-dependent picture of AGNs described above, and our in¬ 
ability to directly resolve the BLRs of even nearby AGNs with current technol¬ 
ogy, it is necessary to constrain the black hole mass scaling factor through some 
means. This is typically accomplished in an indirect way by comparing the re¬ 
lationship between black hole mass and host-galaxy bulge stellar velocity disper¬ 
sion, the Mbh — <7* relationship, for nearby quiescent galaxies with stellar- and gas- 
dynamical modeling-based black hole masses to the relationship for AGNs with as- 
yet-unscaled reverberation masses, i.e., virial products. The average multiplicative 
factor that must be applied to the AGN virial products to bring the AGN relation¬ 
ship into agreement with that_of the quescent galaxies is found to be ~ 4 — 5 (most 
recently 4.3 ±1.1, Grier et ah 2013 a, when all AGNs are treated equally and local 
high-luminosity AGNs are included to extend the range of the Mbh ~ <7* relation¬ 
ship). 

Certainly, there are many assumptions involved in determining this average scal¬ 
ing factor, and it is important to keep in mind when applying the scaling factor that 
it is a population average and therefore likely to be uncertain by a factor of 2 — 3 
for any individual AGN. Nonetheless, several lines of independent evidence point 
to this value of the mass scaling factor being in the right ballpark. 

The first is to simply assume that the main contribution to the scale factor is the 
inclination of the system, or the sin / term in the gas velocity. A scale factor of 4.3 
would imply a typical broad-lined AGN inclination of ~ 29°, which is well in line 
with expectations based on our current understanding of AGN structure. Further¬ 
more, this average inclination agrees well with the geometric inclinations derived 
for the exte nded and resolved narrow-line region structure of several nearby AGNs 


(Fische r et a h. 1201 3). 


A more rigorous test is to compare the black hole masses derived from reverber¬ 
ation mapping with those derived from stellar- or gas-dynamical modeling in the 
same objects. As previously discussed, most AGNs are too far away for dynamical 
modeling techniques to be applied, but a few very nearby AGNs can be examined 
in this way. Direct comparisons have been carried out for two AG Ns to date — 


in tnis way. Direct comparisons nave been c arried out tor two Au [ ns to date — 
NG C 3227 (Davies et all 120061 : iDennev et all l2010h and NGC 4151 (Onken et al 


2014; (see Figure^. So far, the resultant black hole masses agree remarkably well 


between such disparate measurement techniques, each derived from independent 
observations and each with their own independent set of assumptions and biases. 
While the agreement is reassuring for these two objects, a sample of twojs hardly 
definit ive. Two more A GNs with reverber ation masses — NGC 6814 (Bentz et all 
2009bl) and NGC 5273 ( Bentz et al. . 2014 ) — are in various stages of the dynamical 


modeling process, and a handful of other AGNs are being targeted for reverberation¬ 
mapping with the hope of dynamical modeling to follow. 
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Fig. 3 Comparison of black 
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One current complication to this test, however, is the effect of bars on the central 


stella r dynamics, and theref ore the derived black hole masses (e.g.. iBrown et al 


20131 Hartmann et~aD 2014 ). Even for the relatively face-on galaxy NGC4151 


Onken et al. ( 2014 ) found a significant bias was induced in the best-fit black hole 


mass from the weak galaxy scale bar. NGC5273 will be an especially interesting 
case for testing black hole masses from reverberation mapping versus stellar dynam¬ 
ical modeling given its unbarred SO morphological type. Furthermore, in the next 
several years, it is likely that JWST will allow some advancement in the n umbe rs of 


AGNs that can provide direct mass comparisons across techniques. Giiltekin et al 


( 2009 ) argue that dynamical modeling can still place strong constraints on black 
hole mass even if the radius of influence of the black hole is not strictly resolved 
in the observations. JWST will provide a comparable spatial resolution for studies 
of host-galaxy stellar absorption features to that which is currently achieved with 
ground-based observatories and adaptive optics, but the advantages of JWST in¬ 
clude a stable PSF, a significantly higher Strehl ratio, and very low backgrounds, all 
of which are important for deriving tight dynamical constraints on the black hole 
mass. 

Finally, as we discuss in Section ITJl it is possible to directly constrain the 
black hole mass, without needing to resort to the use of a scaling factor, from the 
reverberation-mapping data itself. Data quality concerns have generally not allowed 
this goal to be met in the past, but recent progress is encouraging, and the resul¬ 
tant mass constraints generally agree with our expectations based on the arguments 
above. 
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Fig. 4 The relationship between the H/3 ti me delay and the s pecific luminosity of the AGN at 
5100 A, the Rblr — ^-agn relationship. From lBentz et al.ll2013l) . 


3.2 Black Hole Scaling Relationships 


One of the most useful scaling relationships to arise from the compendium of 
reverberation-mapping measurements is the Rblr — A\gn relationship — the rela¬ 
tionship between the time delay, or average radius of emission, for a specific emis¬ 
sion line and the luminosity of the central AGN at some particular wavelength. 
This particular scaling relationship was expected from simple photoionization argu¬ 
ments and looked for in the early days of reverberation mapping when the number 
of successful monitorin g campaig ns and reverberation measurements was still very 
small ( Koratkar & Gaskellill99 1). The addition of several high-luminosity local PG 
quasars to the reverberation sample led_to the first well-defined functional form of 
the Ablr - Lack relationship dKaspi et all 12000b . Figure® shows the most current 
calibration of the /?blr — Tagn relationship between H/j and the AGN continuum lu¬ 
minosity at rest-frame 5100 A, where all luminosities have been carefully corrected 
for the contamination of host-galaxy starlight using high-resolution HST images and 
two-dimensional decomp ositio ns of the images to separate the galaxy starlight from 


the AGN emission ( Bentz et ah . 20131) . 


Rather than carrying out a long-term monitoring campaign for any AGN of in¬ 
terest, the Rblr — T \gn relationship allows a single spectrum to provide an estimate 
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of the black hole mass through two simple measurements: the width of the broad 
emission line as a proxy for the gas velocity, and the continuum luminosity of the 
AGN as a proxy for the time delay expected in the emission line. This handy short¬ 
cut provides a means for taking large spectr oscopic surveys and producing catalogs 
full of black hole mass estimates (e.g., Shen et al. 20111) . Of course, as one might 
expect, the devil is in the details. 

Currently, H/i is the only emission line for which a well-calibrated 7 ?blr — Lagn 
relationship exists. Unfortunately, H/j shifts out of the observed-frame optical band¬ 
pass at redshifts of only z = 0.7 — 0.8. The Mg II and C IV emission lines in the 
rest-frame ultraviolet are therefore more accessible for most quasar surveys con¬ 
ducted from the ground. But at this time , only a handful of measurements of C IV 
reverberation time delays exist (see Kaspi et ahll2007 for a first attempt to constrain 
a C IV /?blr — Tagn relationship), and even fewer reverberation measurements ex¬ 
ist for Mg II. Quasar black hole masses therefore require bootstrapping the esti- 
mates into the UV using the H j3 7 ?blr — Tagn relationship as the cornerstone (e.g., 
Vestergaard & Peterson 20061) . Furthermore, the exact prescription for turning two 
simple spectral meas urements into an unbiased mass estimate is still highly debat¬ 
able (see, e.g., Denney et al. 20091 2013 for discussion of several of the specific 
details that can cause biases). 

Interestingly, the small scatter in the /?blr — Tagn relationship has led to the pro¬ 
posal that it may be used to turn AGNs into standardizable cand les for investigating 
cosmological expansion ( Watson et al. . 2011 : King et all 2014t) . Quasars are easily 
observed well beyond z»l, where Type la supernovae become rare and difficult to 
find and where the differences in cosmological models are more apparent. One of 
the largest practical difficulties in turning this idea into reality, however, is again the 
long time delays involved in monitoring high redshift quasars and the necessity of 
high signal-to-noise spectroscopy over the course of such a monitoring campaign. 

Other black hole scaling relationships include the aforementioned AGN Mrh — 


g k rela tionship ( Onken et al.L l2004: Gra ham et al.Ll201 ltlPark et al.L l2012: Gri er et al 


2013al) and the relationship between the AGN blac k hole mass and the host- 
galaxy bulge luminosity, the M^h - /-bulge relationship ( Wandel . 2002!: Bentz et al 


2009a : Bentz , 2011). While the more commonly-used forms of these relationships 
tend to be t hose derived for black holes with dynamically-modeled black hole 
masses (e.g., Magomanetalil998 ; Ferrarese & Merrittl2000t Gebhardt et al . 200(1 


Gultekin et al.ll2009HMcConnell & Mall2013l:lKormendv & Holl2013l) . the AGN re- 

lationships provide a useful counterpoint given the differences between the two 
samples. In particular, the AGN reverberation sample has a large percentage of late- 
type, disk-dominated galaxies, whereas the quiescent galaxy sample with dynamical 
black hole masses is comprised mainly of early-type galaxies. It is not at all clear 
that galaxies of different morphological t ypes s hould all follow the same scaling 
relationships (cf. the recent review by Kormendv & Hoi 20 13). Furthermore, active 
g alaxies may not follow the same scaling relationships as quiescent galaxies (e.g., 
WandelH~1999l) . but most studies of galaxy and black hole co-evolution at cosmolog¬ 
ical distances are necessarily limited to active galaxies. 
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Even despite these differences between the samples, it appears that our current 
constraints on the general forms of the Mbh — Tbuige and Mbh — <7* relationships 
are consistent for AGNs with reverberation masses and for quiescent galaxies with 
dynamical black hole masses. Several ongoing studies aim to refine and more ac¬ 
curately constrain these relationships, and part of this effort is devoted to tackling 
the key observational uncertainties that remain — such as determining accurate dis¬ 
tances to the AGN host galaxies, and replacing stellar velocity dispersion measure¬ 
ments from long-slit spectra with those obtained from integral field spectroscopy of 
the host galaxies. Here again, barred galaxies (and galaxies with “disky” bulges) are 
a source of c onfusion. While they are seen to be outliers in the quiesce nt Mbh — <7* 
relationship (Hu, 2008 : Graham et all 2011 : Kormendv & Hoi 2013|h there is no 
such offset seen in the AGN Mbh — (7* relationship ( Grier et all 2013ah unless it is 
artificially inserted by scaling the black hole masses in those galaxies by a different 
value (Ho & Kim, 20141) . Comparisons between the different assumptions and bi¬ 
ases in the AGN Mbh — <7* relationship versus the quiescent galaxy relationship will 
therefore help to shed light on the underlying causes for such puzzles. And while 
it has been the focus of less intense study in the last decade or so, an accurately 
calibrated Mbh — Lbuige relationship will be especially necessary for constraining 
galaxy evolution through upcoming deep all-sky photometric surveys with no spec¬ 
troscopic component, like LSST. 


3.3 BLR Geometry and Kinematics 

Most of the progress in reverberation mapping has focused on the lowest-order 
measurement that can be made, namely the average time delay of an emission 
line, because this is what is required to make a dynamical mass measurement. 
Historically, it was also the only measurement that could be recovered from the 
marginally-sampled light curves from early reverberation campaigns. However, 
there is much more information encoded in densely-sampled light curves. In par¬ 
ticular, the emission-line light curve is a convolution of the continuum variations 
and the extended response of the BLR gas at different line-of-sight velocities and 
light travel times relative to the observer. Resolving the time delays as a function of 
velocity across an emission-line profile can therefore give constraints on the detailed 
geometry and kinematics of the BLR gas. In Ligure[5] we show three examples of 
emission-line response given one fairly simple model for the BLR and three dif¬ 
ferent possible kinematics of the gas: rotation, infall, and outflow. The differences 
between the three are apparent in a full deconvolution of the emission-line response 
(Ligure[5] shaded regions), or a first-order analysis in which the mean time delays 
are computed for velocity bins across the line profile (Ligure[5] error bars). 

While the wealth of information that is potentially available from reverberation¬ 
mapping datasets has been understood for quite some time, the practical difficul¬ 
ties involved in deconvolving a faint signal from sparsely- and irregularly-sampled 
noisy data have limited much progress in this area. Notable early attempts include 
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Wanders et al. ( 19951) . Done & Krolik (1996), and Ulrich & Horne ( 19961) . but re¬ 
covered maps of velocity-resolved results, reminiscent of those shown in Figure[5] 
were ambiguous at best. In the last few years, however, reverberation programs have 
enjoyed much more suc cess given the careful experimental setup and the resultant 
high quality of the data dBentz et al. 2009b . 2010a : Denney et al. 2010t Grier et al 


20lilDe Rosa et al.ll2015lf 

With such high quality data now in hand, there are two general methods for ex¬ 
tracting the reverberation signal. Deconvolution techniques can be used to produce 
a model-independent, but potentially difficult-to-interpret, velocity delay map like 
the models displayed in Figure[5] Direct modeling of the spectroscopic data, on the 
other hand, produces easily-interpreted constraints on different possible physical 
models, but is by definition model-dependent and human imagination-limited. 

The most widely employed deconvolution algorithm to date is the MEMECHO 
code ( Horne et al .1 199 ltlHornel 1994 . which uses maximum entropy balanced by 
to find the simplest possible solution that fits the data. Each pixel in the AGN spec¬ 
trum can be treated as a separate velocity bin for which an entire light curve exists 
throughout the monitoring campaign. The code solves for the time delay response 
function in all the individual light curves for each of the pixels in the spectrum. 
From the many reponse functions, a map of time delay as a function of velocity is 
reconstructed. 





Fig. 5 The expected emission-line response for a toy model BLR with three different possible 
kinematics: (a) rotation, (b) infall, and (c) outflow. For simplicity, the geometry is kept the same 
for all three cases — the line emission is restricted to a bicone with a semi-opening angle of 30° 
and the model is inclined at 20° so that the observer is inside the beam. The radiation structure 
within the BLR clouds is set so that the emission is enhanced for clouds at smaller radii, and 
the line emission is partially anisotropic, such that the emission is enhanced in the direction of 
the illuminating source. The gray-scale images show the full two-dimensional structure in time lag 
versus line-of-sight velocity, while the vertical red error bars show the weighted mean and standard 
deviation of the time lag within discrete velocity bins that are represented by the horizontal red error 
bars. The overall shape is different for each of the three models: a symmetric structure around zero 
velocity for circular Keplerian orbits, longer lags in the blues hifted e mission for infall, and longer 
lags in the redshifted emission for outflow. From lBentz et all li2009b f. 





























































AGN Reverberation Mapping 


13 



4500 5000 5500 6000 6500 


X (A) 

Fig. 6 The deconvolved emission-line response as a f unctio n of veloci ty for the broad optical 
recombination lines in the spectrum of Arp 151. From lBentz et al.N2010ah . 


Figure[6] shows an example of the deconvolved response of the broad optical 


recombination lines in the spectrum of the Seyfert galaxy Arp 151 dBentz et al 


2010a ). The differences in expected mean time delays for the lines is immediately 


apparent, and also interesting is the strong prompt response in the red wings of the 
Balmer lines compared to the lack of prompt response in the blue wings of those 
same emission lines. An asymmetric response such as this could be produced by 
either rotating gas with enhanced emission in one location (such as a hot spot or 
warp), or by infalling gas, or some combination of these simplistic models. Similar 
asymmetries are also seen in the deconvolved responses of a handful of additional 
AGNs dGrier et allfeoi 3bl) . 


Direct modeling, on the other hand, tests the data against specific geometric and 
kinematic models to constrain the family of models that best represent the obser¬ 
vations. For the Arp 151 dataset above, direct modeling results prefer a thick disk 
BLR geometry, inclined at ~ 25° to the observer’s line of sight, and the kin emat 


ics are dominated by inflow with some contribution from rotation ( iPancoast et al 
2014 in general ag reement with the interpretation of the deconvolution results. 


Pancoast et al . (2014) also find similar results for a few additional AGNs. 


Furthermore, direct modeling is able to determine the individual scaling factor 
that would need to be applied to the reverberation-based black hole mass as de¬ 
scribed in Section l3TTI above. For the handful of AGNs with successful dynamical 
models, we can see that the scale factor indeed varies from object to object, as ex¬ 
pected for a population of objects with random inclinations within some limited 
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Fig. 7 The individual mass 2 0 

scale factor for 5 AGNs 
with direct modeling of their 
velocity-resolved responses. j 5 

The average value for this 
small sample is / as 4.8 (de¬ 
noted by the horizontal blue "^ 1.0 

dashed line), in good agree- 

ment with values determined m 

O 

by comparison of the active — 0.5 

galaxy and quiescent galaxy 
AfgH — a* relationships (de¬ 
noted by the o ther horizonta l 0.0 

li nes). From IPancoast et ahl 
(20 |li. 
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range (see Figure[7}. The average of these individual scaling factors also agrees 
quite well with the population average derived above through use of the Mbh — O* 
relationship (see Section ITTV 


4 Looking Ahead 


The future is looking quite busy for applications of reverberation mapping. In addi¬ 
tion to the many ongoing areas of study summarized in the previous pages, several 
recentor ongoing programs, such as the multi-object SDSS reverberation program 


(She n et al.L 120151) and the massive H ST plus ground-based reverberation program 


for NGC5548 (De Rosa et all 2015), are just starting to report results that are sure 


to lead to new insights and new puzzles in AGN physics. Additionally, the upcom¬ 
ing OzDES program (King et all 201511 will help push to higher redshifts, providing 
a more stable anchor for black hole mass estimates of high-z quasars. 

The flurry of recent activity in velocity-resolved reverberation mapping is un¬ 
likely to abate any time soon, and here we may hope to unlock many of the secrets 
surrounding AGN feeding and feedback. New codes to deal with velocity-resolved 
reverberation mapping data are currently being developed ( Skielboe et al.l2015 , An¬ 
derson et al. in prep), and new features are being added to currently-existing codes 
(Pancoast, private communication). 

With UV astronomy currently dependent on the contin uati on of //.ST, it is cer¬ 


tainly conceivable to think that the Kronos spacecraft (Peterson et al., 20031) . or a 


similar instrument, may again make an appearance in proposal form. It is an excit¬ 
ing time for supermassive black holes, reverberation mapping, and AGN physics! 
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